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Topological phases which host Majorana fermions cannot be identified via local order parameters. 
We give simple non-local order parameters to distinguish quasi- ID topological superconductors of 
spinless fermions, for any interacting model in the absence of time reversal symmetry. These string 
or brane-type order parameters are amenable to measurements via quantum gas microscopy in cold 
atom systems, which could serve as an alternative route towards Majorana fermion detection that 
involves bulk rather than edge degrees of freedom. We also discuss our attempts at accessing 2D 
topological superconductors via the quasi- ID limit of coupling N identical chains with Zjv translation 
symmetry. We classify the resulting symmetric topological phases and discuss general rules for 
constructing non-local order parameters that distinguish them. The phases include quasi-lD analogs 
of (i) the p + ip topological superconductor, which can be distinguished up to the 2D Chern number 
mod 2, and (ii) the 2D weak topological superconductor. The non-local order parameters for some 
of these phases simply involve a product of the string order parameters for the individual chains. 
Finally, we sketch a physical picture of the topological phases as a condensate of Ising charged 
domain walls, which motivates the form of the non-local order parameter. 



Quantum phases with emergent Majorana fermion ex- 
citations have received much attention in the past sev- 
eral years ^ with recent experimental evidence suggest- 
ing they have been realized in solid-state systems^ Ma- 
jorana fermions are known to appear at the boundaries 
and topological defects of exotic ID 3 - and 2D topolog- 
ical superconductors.^ Cold atom realizations of such 
phases would serve as a new platform for studying and 
manipulating Majorana fermions. 

At the same time, a general framework for classify- 
ing quantum phases continues to be developed. Im- 
portant achievements include the classification of free 
fermion systems^^ an understanding of how inter- 
actions can modify these results in certain symme- 
try classespT— and general methods for many of the 
symmetry-protected bosonic or fermionic systems with 
interactions J^r— In ID, where matrix-product states 
provide a framework for describing ground states, the 
gapped symmetry-protected topological phases have 
been completely classified! 11 ' 13 A natural resulting ques- 
tion is how such phases may be distinguished via accessi- 
ble quantities, in any dimension. Fully symmetric phases 
are distinct but have no broken symmetries and hence 
are immune to a local order parameter description. 

The insight obtained from classification enables de- 
sign of non-local order parameters that extract the defin- 
ing quantities, associated with the cohomology group of 
the symmetry group, which characterize a symmetry- 
protected phase. This problem has been largely ad- 
dressed in ID for boson o 17 ' 18 but remains open in higher 
dimensions. These ID bosonic non-local order parame- 
ters (OPs) are robust in that they are strictly symme- 
try rather than wavefunction dependent; their extension 
to fermions to our knowledge has not been discussed. 



In contrast, conventional AKLT string order 
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vant for the Z2 x Z2 spin rotation symmetry protected 
Haldane phase of ID bosonic systems, measures extra- 



neous aspects of the wavefunction in addition to topo- 
logical information. Nonetheless, AKLT string order has 
proven useful in some contexts. Analogously, we address 
whether it may be possible to construct similar string or 
brane-type^ non-local order parameters for other sys- 
tems. Attempts at extending string order to brane or- 
der for systems with dimensionality greater than strictly 
ID have been discussed before ) 22 ' 23 in particular in the 
context of the Haldane phased In contrast, our basic 
building block is the Majorana chain. 

We first consider quasi- ID spinless fermion topological 
superconductors with interactions but no symmetries. It 
is known that AKLT string order can be derived via a 
non-local mapping which transforms topological proper- 
ties of the Haldane phase into broken symmetries , 20 1 2 1 1 25 
Likewise, the non-local Jordan- Wigner mapping can 
transform certain fermionic topological properties to 
bosonic symmetry breaking. This yields non-local OPs 
distinguishing the topological and trivial phases of our 
system; the bulk of these order parameters measures 
fermion parity, while their terminating operators may ei- 
ther be fermionic or bosonic. 

Our order parameters are relevant for cold atom exper- 
iments, which have seen recent breakthroughs with the 
development of the quantum gas microscope ) 26 ' 27 as well 
as subsequent measurements of non-local order ^ by the 
groups of M. Greiner and I. Bloch, respectively. In these 
measurements, the parity of the particles is simultane- 
ously measured on all sites of the lattice, which consti- 
tutes a non-local probe of the many-body system that is 
particularly well suited to identifying topological phases. 
In contrast, most other probes measure local properties, 
such as correlations between a pair of local operators. 
While these local probes are well suited to identifying 
phases with a local order parameter, they are blind to the 
subtle correlation in the bulk of topological phases. Mea- 
suring topological aspects of free fermion band structures 
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in cold atom systems have been discussed^— - here, 
however, we will be concerned with generic interacting 
systems, in particular topological superconductors. As 
an example of a quantity measurable with current exper- 
imental techniques, we describe a system of two identical 
chains for which a topological phase can be detected via 
fermion parity only. 

We then add Zjv translation symmetry to N identical 
chains with interactions and describe the protected sym- 
metric phases, which can scale to interesting 2D phases. 
The classification appears to distinguish certain topolog- 
ical indices in the case of free fermions. We justify simple 
symmetry transformation rules which the edge operators 
of a non-local OP should satisfy in order to uniquely dis- 
tinguish the symmetric phases. As further illustration, 
we discuss a condensate of Ising defects. 

Our results can be supported by either working in 
bosonic or fcrmionic variables; in the main text we mainly 
take the bosonic point of view and discuss the fermionic 
description, also sketching a derivation of "selection 
rules" for non-local OPs, in the appendices. Through- 
out the paper, we use the term non-local OP to include 
string and brane-type order. 



I. INTERACTING SPINLESS FERMION 
TOPOLOGICAL SUPERCONDUCTORS 

Example: Single Majorana Chain. To illustrate 
the general form of the order parameters, we first consider 
Kitaev's spinless p-wave topological superconductor on 
an open chain with Hamiltonian: 

H Kit = y^(-ta|aj + i + \A\a,ia i+1 + h.c.) - ^(ajaj - 

(la) 

= \ XX H + l A Dx*+i + (* + |A|)XiXi+i - MXiXi} 

(lb) 

with site fermions ai = \{xi + *Xi) an d Majorana opera- 
tors Xii Xi^ The phase of the superconducting order pa- 
rameter A = \A\e t6 has been gauged away. For |^| < 1, 
there are gapped topological phases if |A| ^ and a 
gapless normal phase if |A| = 0; if |£| > 1 for any 
A|, the phase is gapped and trivial. Let the Jordan- 
Wigncr mapping be erj = e* 7m -' , crj = n»<j e ™ ni Xj, a j = 
~ T\i<j e%lxni Xj- The Majorana chain maps onto an XY- 
type spin model in a transverse magnetic field. Fermion 
parity corresponds to a Z 2 spin symmetry . erj . This 
symmetry is broken or unbroken, respectively, in the 
spin system when the corresponding fcrmionic system is 
gapped and in a topological or trivial phase. This is in 
fact a more general correspondence between ID fermionic 
and bosonic systems. 

Consider an Ising limit (e.g. set t = |A|) of the spin 
Hamiltonian obtained from the Kitacv model: H SV i„ = 



Ei + 9°*) with J = |A| and g = A 
two-point spin correlation function is finite in the spin 
ordered phase and vanishes in the disordered phase. It 
maps to a string OP which distinguishes the topological 
phase from the trivial in this limit: 

fe-l 

«^fe) = (HXi) II e "" 3 ») ( 2 ) 
j=i+i 

As fermion parity cannot be broken, the fermion ground 
states in the topological phase correspond to linear com- 
binations of the 1i breaking spin ground states. Two- 
point correlations are, however, insensitive to the choice 
of linear combination. 

The Kitacv model Eq. [T]has additional symmetries, for 
instance time reversal x ~ * XiX ~X'i this reduces the 
possible two-point spin correlations which can be chosen. 
a z or cr v correlations are finite in the broken symmetry 
regimes of t > or t < 0, respectively, so that we have 
either (XijXfc) string termination operators, as in Eq. [2j 
or (xi,Xfe). For fermion models with strictly no other 
symmetries, these constraints will not occur. 

On the other hand, the self-duality of the quan- 
tum transverse Ising model under mapping t x L = 

f°f+i' T i+i = n^i "! to domain wall variables on 

bonds yields H du ai = Y,i(-\^\ T i + L + \^ T i^L T i + i) in 
the thermodynamic limit. A two-point correlation in the 
r variables distinguishes the two phases^ In the t — |A| 
limit, this yields a fermion string OP which is finite in 
the trivial phase and vanishes in the topological: 

k k 

General Form. We used the Kitaev model, and in 
particular, an Ising limit, to illustrate a more general 
correspondence which holds for quasi-lD topological su- 
perconductors of spinless fermions with interactions and 
no symmetries. These fermionic phases were classified 
in Refii by considering the bosonic phases protected by 
a global bosonic Z 2 symmetry corresponding to fermion 
parity. There are only two gapped phases, which have 
the Z 2 symmetry unbroken or broken; via the Jordan- 
Wigncr mapping, they correspond to fcrmionic symmet- 
ric phases that are, respectively, trivial or topological, 
with boundary Majorana zero modes. The models we 
cite have translation symmetry along their infinite di- 
mension, which multiplies the number of possible phases 
by a factor but we neglect this, focusing only on topo- 
logical distinctions. In this case, there are two distinct 
phases. 

We can distinguish the phases in the bosonic variables 
and map the result to fermions. Any two-point corre- 
lation function (OiO'j) with O/O' local operators which 
are odd under the Z 2 spin operation is generically finite 
as \i — j I —¥ 00 in the the spin ordered phase and vanishes 
in the disordered phase. Hence, the two-point function 
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maps to a fermionic string OP (for one chain) or brane 
OP (for two or more) whose bulk measures fermion parity 
and which is terminated by fermionic operators. 

The spin disordered (i.e. symmetric) phase is not 
susceptible to local order parameters. Rather, utilizing 
bosonic selection rules proposed ir*i£, we conclude that a 
non-local OP which is finite in this phase should apply 
the local Z2 symmetry over a domain in the bulk, and the 
domain should be terminated by operators which are Z2 
invariant. In Appendix I, we argue that this order param- 
eter vanishes in the ordered phase of spins. Mapped to 
fcrmions, an order parameter which is finite in the trivial 
phase and zero in the topological phase would consist of 
a bulk which measures fermion parity and is terminated 
by local bosonic operators. 

The appearance of fermionic or bosonic terminations 
for a non-local OP is a fermionic selection rule, analo- 
gous to those described in 18 for bosonic systems, which 
distinguishes the two phases. As an alternative to the 
Jordan- Wigner mapping, in Appendix III we try to jus- 
tify fermionic selection rules for non-local OPs from 
fermions directly based on ideas from fermion classifi- 
cation. 

To summarize, order parameters for the two topo- 
logically distinct phases of interacting spinless fermion 
topological superconductors can be constructed with the 
form: 

S top = {0 FL \{e^O FR ) (4a) 
jen 

Striv = (Obl J] ^0 B r) (4b) 

where Oft,/r,Obl/r & re local fermionic or bosonic op- 
erators near the left, right edges of region O. Stop is non- 
zero in the topological phase and vanishes elsewhere; the 
behavior of Striv is reversed. These are the generic val- 
ues, as we now discuss. 

Remark. While the order parameters proposed 
throughout this paper can be used for general interact- 
ing models, their values depend in part on the state, as it 
requires evaluating matrix elements of certain local oper- 
ators. This is no different than tailoring an order param- 
eter for a symmetry breaking theory: certain operators 
may be more "optimal" for detecting the broken symme- 
try because they yield larger values, while specific models 
may have larger symmetry groups, which we can identify 
from the outset. For instance, the symmetry group of 
the quantum Ising model with a z nearest-neighbor cou- 
plings includes time reversal of spins followed by ir rota- 
tion about y, so that (erf) = 0. Our order parameters 
distinguish the phases simply, theoretically and experi- 
mentally, and give an interpretation to the distinct order 
in the phases. Their general form, which is robust to 
interactions, will be given by Eqs. U 

Microscopic Picture. We explain why fermionic 
or bosonic operator ends distinguish the topologically 
distinct fermionic phases. To illustrate, we specialize 





W w KJ WW yJ 


n 


X 


Fu A< hi Pu K, 

S 


X 

1 








00 


OO! •• 


OO 


X X 

1 


XX XX 

1 s 1 

trivial 


X X 



FIG. 1. Single chain at top shows two Majorana fermions 
X, X P er she (circles with same color) with non-trivial pairing 
(boxes). The non-local order parameter Stop for the topolog- 
ical phase measures this pairing (green line). Bottom chain 
shows a phase with trivial Majorana pairings, which is mea- 
sured by Striv 



to string order in the Ising limit of the single Kitaev 
chain. Introduce bond fcrmions dj = + *Xj) 

of re-paired Majoranas, neglecting the non-local fermion 
ajvi, = ^(xi +iX-/v) by working on an infinite chain. This 
basis exactly solves the t = |A|, [i = limit. The topolog- 
ical and trivial phase string OPs ©,(© can be rewritten 

(k >i + iy. 

fc-i fc-i 

Sto P = (-i^) II <'' T " \> * \ [' (5a) 

j=i+l j=i 
k k-1 

Strw = n eiwnj <* (Si-i + n e "" 3 - 4) 

j = i j=i 

(5b) 

Evidently, the fermionic or bosonic nature of the ter- 
minations depends on the basis used. The topological 
ground states at t = | A | , /i = have uniform filling in 
the bulk in the bond fermion basis hi so that St op , which 
measures their parity, is finite. Weak perturbations \i ^ 
drive the system away from this "bond-centered" order- 
ing; they favor on-site Majorana pairings and hence cre- 
ate localized pairs of bond fermion defects with respect to 
the unperturbed state. But because the defects are pairs 
and are localized, they weakly modify the bond fermion 
parity Stop measured in a region. They are more likely to 
fall into the bulk of Stop, in which case they do not mod- 
ify the bond fermion parity measured, rather than cross 
its ends. Hence, the value of Stop remains finite in the 
topological phase. A dual picture holds for site fermions 
a,i deep in the trivial phase, for which site fermions are a 
good basis to use (i.e. the wavefunction is simple). This 
explains the finite value of Striv in the trivial phase. 

On the other hand, such string OPs vanish in the com- 
plementary phases. To understand how this occurs, con- 
sider perturbatively evaluating Stop in the trivial phase 
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|A| Kitacv model. The ground state at the 
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with [i < is the site fermion 



vacuum |0) which is then corrected by the perturbation 



^ = |A|£ ? 



<j l XjXj+i = l A l - Uj)W+i t -j+i 

corrects |0) by creating localized pairs of "defects" rel- 
ative to |0); these become delocalized, or new ones are 
created, with higher orders of perturbation theory. But, 
locality is a strong constraint on physically allowed op- 
erators; V preserves fermion parity not just globally but 
also "locally," in a certain sense. On the other hand, St op 
must connect states which differ in site fermion parity at 
two widely separated points i, k. Such states cannot arise 
through the effects of a local and fermion parity preserv- 
ing perturbation applied to an initial state with uniform 
parity throughout. Hence, St op should asymptotically 
stay zero away from the point Hq. We anticipate that 
this argument holds for the entire phase because of the 
strong constraint of locality on physically allowed fermion 
operators. Striv can also be argued to remain zero in the 
topological phase by a similar, dual argument. 
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II. NON-LOCAL ORDER AND QUANTUM GAS 
MICROSCOPY 



Enabled by advances in single-site resolved imag- 
ing of optical latticeS ) 26 i 27 non-local measurements in 
cold atom systems arc now possible and were re- 
cently demonstrated 2 ^ for string order in bosonic Mott 
insulatorsi 24 i 34 Similarly, 



Striv — 6* 



(6) 



which yields a finite value in the trivial phase, is accessi- 
ble with current experimental techniques. 

We consider how one might measure OPs with more 
complex terminations. For instance, 
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(7) 



directly detects the topological phase by generically 
yielding a finite value. The difficulty with measuring 
non-local OPs such as St op is that they are off-diagonal 
in the site fermion basis imaged in experiments. We sug- 
gest a scheme for measuring a string OP such as S top on 
the interval [f , N) in the bulk of a long Majorana chain. 
The idea is that by evolving the ground state in a con- 
trolled manner, such as with fermion tunneling, we may 
extract the additional information needed to reconstruct 
the string OP (Fig. H]). 

For instance, let a Kitaev chain ground state be = 
Sijfe Pijk\ni)\nf)\nk)- Here, \n{) is a site fermion config- 
uration indexed by i for the inner region sites 2 to N — 1 , 
n^) indexes states for the region outside [f,7V], while 
\rik) is a configuration for the string end sites 1 , TV, with 
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FIG. 2. Top figure shows single Majorana chain (sites are 
circles) and the order parameter Eq. [6] finite in the trivial 
phase, which can be measured with current cold atom tech- 
niques since it only involves fermion parity (highlighted green 
circles). Bottom shows a potential experimental scheme for 
measuring order parameters such as Eq. [7] with complex ter- 
minations. In the bulk (green), fermion parity is still mea- 
sured, but additional measurements for the end sites labeled 
1, N (blue circles) must be made (see text) to extract the 
string order parameter value. 
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value is 
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{|0>, «1 10>, o]v|0), 0^0,^10)}. The measured 



2P l (-Re((3 ijl p ij4 ) + Re{p lj2 fa j3 )) (8) 



where Pi is the parity of configuration i for sites 
[2,7V— 1]. The additional information, beyond ampli- 
tudes |/?r//c|, needed in order to reconstruct the expec- 
tation value is certain relative phases, such as those in 

We imagine consistently starting the system in a fixed 
Majorana chain ground state \ip). A tunneling Hamil- 
tonian Ht which for instance couples only sites 1, TV is 
turned on rapidly, preserving the state. We may consider 
changing the experimental geometry to have the single 
chain folded into two in order to couple I, TV. After dy- 
namic evolution with Ht, the site fermion occupations 
are measured at specified times. Accurate knowledge 
of the coupling Hamiltonian parameters and amplitudes 
|/3yfe|, determined from previous repeated measurements, 
may enable extraction of the necessary relative phases 
and reconstruction of the string OP value. 

The practicality of the suggested scheme for current 
systems would need to be determined. A general chal- 
lenge appears to be the number of measurements needed, 
as a ground state for N sites in the deepest regions of the 
topological phase consists of an exponential in N number 
of states in the site fermion basis, all with equal magni- 
tude weights. Design of a protocol to enable extraction 
of the off-diagonal interference terms, if possible, would 
be interesting. 
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Example: Two Identical Chains. We consider the 
case of two identical Kitaev chains A, B each with param- 
eters (t, \A\,fi) and coupled with an interchain hopping 
t_i_. This high symmetry model is a special case of the 
general N chain system with Zjv translational symmetry 
considered in the subsequent sections, but we emphasize 
it here because of its potential experimental relevance. 

The phases of this system can be easily seen by switch- 
ing to momentum k y = 0, n in the transverse direc- 
tion. The resulting Hamiltonian consists of two decou- 
pled Kitaev models for the k y = 0,7r degrees of freedom 
(DOF) {ao(i)}i U {a„(i)}i with modified chemical poten- 
tials fi± = fi ± 2t±: 

H = Hkh.a + H Kit , B - 2t± ^2(a] A a lB + h.c.) (9a) 

i 

= H Kit fi{n + ) + H Kit ^{v-) (9b) 
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The phases of the system for |A| ^ have two, one, 
or zero Majorana zero modes per edge as the interchain 
coupling t± is increased (the boundaries are the same as 
those in Fig. [4]). The phase with two Majorana zero 
modes per edge is protected by the exchange symmetry. 
To distinguish the phases we need only independently 
test whether the k y = 0, ir DOF are in the topological or 
trivial phases using Kitaev model string OPs. In regimes 
where only one of k y = Q,ir DOF are in the topologi- 
cal phase, we use string termination operators such as 
Xo/x-n- ~ XA ± Xb or those built out of x operators. 
When both k y = 0, tt DOF are in the topological phase, 
the terminating operators are for instance XoXtt ~ XaXb- 
In other words, two copies of the topological phase order 
parameter of the Kitaev model, one each for the k y = 0, n 
momentum DOF, detects the weakly-coupled regime of 
this two chain system, in which each end has two Ma- 
jorana zero modes. This is equivalent to a product of 
topological string OPs for each chain: 



k-i 



j=i+i 
k-l 



(10) 



eu(0*u(i) II e l<nA[j)+nB[j)) XA(k)xB{k)) 

3=t+l 



Taking products of string order parameter works here 
because of the additional protecting symmetry. 

Alternatively, we ask whether it is possible to devise an 
order parameter for the two chain system which involves 
only fermion parity but which detects a non-trivial phase, 
with Majorana zero modes. That this might be possible 
is suggested by the form of Eq. [TU1 which is terminated 
by bosonic operators such as XaXb- In fact, we will later 
give selection rules which the terminating operators of a 
non-local OP should obey in order to uniquely detect a 
symmetric phase among other symmetric phase (see Ta- 
ble UTJ) . For the phase with two Majorana zero modes 
per edge, the terminating operator should be bosonic 



FIG. 3. Geometry of the non-local order parameter Eq. HOI or 
1111 for a system of two identical chains A, B (white circles are 
fermion sites). Bulk (blue) of the order parameter measures 
fermion parity while specially chosen terminating operators 
act on sites of the two chains (orange rectangles) separated 
by a large distance \i — j\. For instance, these terminating 
operators can be taken to be the fermion number difference 
71 a — tib (see text). This non-local order parameter uniquely 
identifies the non-trivial phase with two Majorana zero modes 
per boundary which is protected by the exchange symmetry. 



(even under parity) but odd under exchange. Opera- 
tors with other transformation rules under the symme- 
tries (fermion parity and translation) detect the other 
symmetric phases. 

fc-i 

([n A (i) - n B (i)} J] e -(^«+n B 0-)) [nA{k) _ nB{k)]) 
j=i+i 

(11) 

has the proper transformation rules to detect the phase 
with two Majorana zero modes per edge. This non-local 
OP works, for instance, for a model of two identical 
chains with intrachain pairing and interchain diagonal 
hopping. While it vanishes for the special model Eqs. [§] 
because of its larger symmetry group, we expect that for 
models with no additional symmetries this order param- 
eter will detect the special topological phase. 



III. PHASES WITH ADDED 1 N 
SYMMETRY 



TRANSLATION 



N Chain Systems. Consider a system of N iden- 
tical, interacting topological superconducting chains of 
spinless fcrmions with Zjv translation symmetry trans- 
verse to the infinite chain length. The geometry is that 
of a cylinder with finite circumference N. We first seek to 
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understand the symmetric fermionic phases, which can- 
not be distinguished from each other through local mea- 
surements. One way to identify the gapped quantum 
phases is to classify the corresponding bosonic phases 
using results from the group cohomology approach to 
classification.— Alternatively, the fermionic phases can 
be identified directly by analyzing fermionic symmetry 
operators, as ir>i£. We give a correspondence between 
the two descriptions specific to Z2 x Z^v- We then illus- 
trate with models for the phases. 

Bosonic Classification. ID gapped bosonic phases 
with local interactions are in correspondence with the 
unbroken subgroups G' of symmetry group G and their 
second cohomology group H 2 (G', C).±±>±2 H 2 (G',C) is 
also the group of equivalence classes of projective repre- 
sentations U of G' with factor systems ui € C, meaning 
U(gi)U(g 2 ) = u(gi,g 2 )U(gig 2 ). Projective representa- 
tions U, U' with factor systems uj,uj' are equivalent if 
they differ by redefinitions of the representations, that 
is, if U'(g) = f3(g)U{g) with /3 e C. Certain com- 
plex phases associated with a projective representation 
are invariant under these gauge changes and are there- 
fore characteristic of an equivalence class and of a quan- 
tum phase. These gauge-invariant quantities distinguish 
the symmetry-protected topological phases. The bosonic 
phases we obtain may be fully symmetric or involve sym- 
metry breaking; as we will see, both classes can map to 
fully symmetric fermionic phases. 

We take a winding path in the Jordan- Wigner map- 
ping to transform the fermionic cylinder into a ID infinite 
bosonic chain; each unit cell accounts for one cylindrical 
ring and contains N spin-1/2 degrees of freedom (DOF). 
The 1m fermion translation symmetry for the circum- 
ference, Oij — > mod TV: maps to a Z^r symmetry 
internal to the unit cell, which we note is not translation 
of the spins. 

The possible symmetric bosonic phases depends on 
the parity of N, since H 2 (Jj2 x Zjv,C) = Z2 or Zi for 
N even or odd. To see this, consider a translation- 
invariant matrix product state (Tj , A) with physical DOF 
per site indexed by j. If the state is symmetric under a 
symmetry £(<?) with g £ G, the matrices transform as 
EjvE^jj/IY = e ie »UlTjUg&&. For our Z 2 x Z N sym- 
metry generated by the bosonic versions of fermion par- 
ity P and translation T, the overall phase of the genera- 
tor projective representations Up, Ut can be gauged such 
that Up = U£ = 1, but the phase in U P U T = e^U T U P 
cannot be eliminated and must satisfy e 2 ^ = e Nl< ^ = 1. 
Hence, of the solutions <fi = 0, ir, the latter is forbidden 
for N odd. These gauge-invariant phases are quantized 
and so are preserved under smooth, gap-preserving de- 
formations to the wavefunctions. A product state can 
be represented by Fj and hence Up,Ut all scalars, so 
= describes the trivial phase. <fi = it characterizes a 
topologically non-trivial phase. 

The relevant symmetry breaking bosonic phases are 
those identified with the proper cyclic subgroups gener- 
ated by translation T and also (for N even) the product 



of fermion parity and translation PT when mapped to 
spins; we denote these by G' = (T) and G' = {PT) re- 
spectively. While parity symmetry P is broken in these 
bosonic phases, it must be restored for fermions; these 
two classes retain enough symmetry so that, as a conse- 
quence, all other broken bosonic symmetries in Z2 x Zjy 
are also restored for fermions. That is, the resulting 
fermionic phases are symmetric. 

For a complementary description of the phases, we 
also apply the approach developed ini^ for fermions di- 
rectly. On a certain subspace, fermionic symmetry oper- 
ators acquire an effective form (P, T) consisting of two 
"fractional" pieces; again, it is the commutations of these 
pieces which identify the quantum phases. We elaborate 
on this in Appendix II for our case and mention the re- 
sult here. Two angles fi, fi' , defined such that if symmetry 
operators acquire effective forms P ~ P^Pp,T ~ T^Tp 
then P L P R = e^P R P L and T L T R = e l ^T R T L , are suffi- 
cient to characterize the quantum phases. 

The correspondence between the bosonic and fermionic 
classifications is given in Table HI In general, fractional- 
ization of fermionic P into fermionic pieces means it is 
broken in bosonic variables^ This is not true for other 
fermionic symmetries (such as Zjv translation) whose be- 
havior in bosonic variables depends in part on those of 
parity (see Appendix II). 



Bosonic 


Fermionic 


G.S.D. 


Physical 


description 


description 




example 


1. Trivial 


(0,0) 


1 


Trivial 


symmetric G' — G 








2. Non-trivial 


(0,0 


4 


Weak Top. 


symmetric G' = G 






Supercond. 


3. Symmetry breaking 


(tt.O) 


2 


Strong Top. 


G' = (T) 






Supercond. 


4. Symmetry breaking 


(n,n) 


2 


Strong Top. 


G' = (PT) 






Supercond. 



TABLE I. Quasi-ID fermionic symmetric phases, henceforth 
labeled 1-4, for N even with symmetry group G—Z2 x Zjv. 
Ground state degeneracy (G.S.D.) listed is for the generic case 
with open boundaries. For N odd only Classes 1, 3 exist. 



Representative Models for Symmetric 
Fermionic Phases. We give models for which 
the symmetry operators take simple effective forms, 
thereby allowing matching to the correct phase. 

The fermionic operators P and T are: 

P = TT e l ™"y W T = TT e tky7lk y (i) (12) 

where i, k y respectively index lattice sites along the cylin- 
der length (x) and momentum around the circumference 
(y) . nk y (i) measures the occupation of the mode au y (i) = 

^= 2j=i £ %ky ^&i,j- A decomposition into Majorana op- 
erators Xk y («): Xk v (i) is a ky (i) = \ [x ky (i) + iXk v (»)] ■ 
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We can view the set of operators {afe„(«)}i for fixed k y 
as DOF for a single Majorana wire. To construct a model 
for Class 3 (symmetry breaking bosonic phase with G' = 
(T)), we consider fixing the ground state occupations of 
all the k y ^ DOF (i.e. put all such k y ^ chains 
into the trivial phase) so that T will act as a scalar in 
the ground state subspace. By treating the remaining 
k y = DOF as a topological Majorana wire, P splits 
into Pl,Pr components which are fermionic. Therefore 

H 3 = 5Z«Xo(j)xo(« + l)+ J2 (*)**,(*) (13) 

serves as a representative for this class. The symmetry ef- 
fective forms on a subsystem spanning sites i E [1, L] of a 
closed system are P = ixo{^)xo{L),T = 1. On the cylin- 
der, systems in Class 3 generically have a ground state de- 
generacy which can be determined from symmetry break- 
ing for bosons, namely \G/G'\ = 2 with G' = (T). 

To construct a model for Class 4 (for even N) , we sim- 
ply switch the treatments of the k y = 0, tt DOF. 

H4 = ^iXn{i)x*(i+l)+ iXkyiijXkyii) (14) 

serves as a representative. Class 4 systems on the cylin- 
der also have two-fold ground state degeneracy generi- 
cally since \G/(PT)\ = 2. 

We can construct a model for Class 2 (for even N) by 
placing both the k y = and k y = tt DOF into the topo- 
logical phase and fixing the site occupations of the re- 
maining k y 7^ 0, 7r DOF. A representative model is there- 
fore 

H2 = X^X- + 1) + *Xo(*)xo(* + 1)) + 

V ■ n- n (15) 

For a subsystem spanning sites i G [i-,L] of a closed 
system, the effective forms on the subsystem are P = 
[iXo(l)xAl)}[iXo(L)x*(L)},f = i X M)Xn(L), so that 
(//,//) = (0,7r). This particular model has four-fold 
ground state degeneracy, but it is physically plausible 
that Majorana zero modes can be gapped out in pairs on 
the edges without the system undergoing a topological 
transition. We expect that four is the minimal ground 
state degeneracy of this phase. 

Physical Models. The previous models become in- 
creasingly non-local for large N; we connect them to 
phases of a local non-interacting model. A nearest neigh- 
bor interchain hopping t±_ and, for N > 2, a nearest 
neighbor interchain pairing Aj_ are allowed by transla- 
tional invariance. We consider the simple lattice p + ip 
topological superconductor (TSC) studied ir>^. The 



Hamiltonian is (N > 2, considering N even): 

j=N 

H = (-taj tj cn+i,j + |A|oi,jOi + ij + h.c.) 

i,j=l 

-Hiriij - -) + (Aj_ dijCtij+i - t±a\jai 1 j+i + h.c.) 

(16) 

with i,j indexing sites along the cylinder length and 
circumference, respectively. We take fixed parameters 
\Im(A±_)\ > 0,i?e(A ± ) = and |A| > 0. In this case, 
there are transitions between phases including quasi-lD 
versions of the 2D weak and strong TSCs as — , — are 

varied.— 

The phases of the system can be seen by rewriting the 
Hamiltonian (fTB|) as 

H = Hku,o(m+) + H K itA»-) + h 'ku (17) 

where Hku,o, Hku^ are the Kitacv Hamiltonians (fTJ) 
with (xo,Xo) an d (x?r, A>) Majorana DOF, respectively, 
and, as before, fj,± = [i ± 2ij_. To analyze the remaining 
piece H' Kit containing all k y ^ 0,ir DOF, we transform 
the Majorana basis by recombining the four Majoranas 
for each fc = \k y \ ^ 0,tt into r) ko /5 ko (i) = ^(x&oW ± 

X-feo(O) and Vk /8k (i) = ^(X-feoW ± Xfeo (*'))■ Thc 
remaining Hamiltonian can be viewed as a collection of 
two Kitaev chains for each kg 6 (0, tt) with Majorana 
DOF (77, 77) and (S, 5) and modified chemical potential 
fik = + 2ij_ cos(fco). There is an "interchain" coupling 
in this basis which is proportional to 7m(Aj_) and which 
gaps out the Majorana zero modes of each chain: 

feoe(0,7r) 

- Um(A ± )sin(k )(r]k Sk +Sk f]k )(i) 

fcoe(0,7r),i 

(18) 

Consequently, the phases of the system are determined 
by whether the independent k y = 0, tt DOF are topolog- 
ically non-trivial. The k y = DOF are in the topo- 
logical phase when the associated chemical potential is 
sufficiently weak \t± + ^ < \t\. Here the ID Z2 invariant 
v ky =o = 1<^ Likewise, the k y = ir DOF are topological 
when Vk v =TT = 1 for \t± — ^ < \t\. Hence, for weak |^-| 
and |-| > i the system has two Majorana zero modes 
per edge, at k y = 0, tt. This phase scales to a 2D weak 
TSC as N ->■ 00 in Eq. ([T6]). as in this regime the 2D Z 
invariant (Chern number) is v = 0. For intermediate val- 
ues |fj_| ~ \t\ and weak chemical potential the system has 
a single Majorana zero mode per edge (either at k y = 
or tt) and will scale to a 2D strong TSC as N — > 00 since 
the 2D Z invariant \u\ = 1. When interchain hopping ij_ 
dominates over intrachain hopping | — | > |^| + 5, the 
system is a weak TSC in the x-direction. Fig. U gives a 
phase diagram. 
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Since the three sets of DOF k y = 0, k y = ir, and 
k y 7^ 0, 7r decouple, we can tune each into topological or 
trivial phases independently while maintaining transla- 
tional invariance. For instance, if the k y = DOF form 
a non-trivial state, treating this as a single chain with 
no other symmetries, we can find a path connecting to 
the model J^. *Xo(*)Xo(* + 1) which preserves translation 
symmetry since it only involves k y = operators. Hence, 
the quasi-lD phases of Eq. [16j namely the 2D weak TSC 
associated with y-direction layering and the two strong 
TSCs with {v ky=0 ,v ky= ,) = (1,0) or (0,1), would fall 
into Classes 2, 3, and 4, respectively, of our classifica- 
tion. It appears that for free fermions our classification 
identifies the ID Z2 invariants Vk y =o and Vk y =-n and con- 
sequently v mod 2 since Vk y =o + Vk y =ir = v mod 2^ 

The weak TSC with Vk x i= for k x = or 7T, which 
is associated with layering in the x-direction, appears as 
a trivial phase. This phase results with strong t±, but 
in the momentum k y basis this coupling is an effective 
chemical potential, \i± = [i ± 2t±, which favors on-site 
pairing of y-momentum Majoranas, driving the k y = 0, n 
Majorana chains away from non-trivial pairing. It is nat- 
ural that that our classification is unable to detect the 
topological index associated with x-translation symme- 
try along the cylinder length, as only y-translation has 
been included. 



IV. NON-LOCAL ORDER PARAMETERS FOR 
THE SYMMETRIC PHASES 

We construct non-local order parameters to distinguish 
the Z2 x Zjy protected symmetric phases of the previous 
section from each other; we consider even N , which en- 
compass the results for odd N. The construction is based 
on general distinctions made apparent by the classifica- 
tions; working in the ID thermodynamic limit, we use 
the bosonic point of view for illustration. 

The bosonic description has revealed hidden struc- 
ture (breaking of certain symmetry operators) which can 
be used, along with recently derived selection rules for 
bosonic non-local OPs^ to identify bosonic operators 
which distinguish the four phases and then map back to 
fermions. 

In the infinite bosonic chain, there is a natural unit cell 
which makes the Zjy translation symmetry on-site. To 
distinguish symmetric bosonic phases, we chose a symme- 
try to apply over many unit cells of the chain (a string), 
and we terminate the domain with operators obeying 
proper symmetry transformation rules. Mapped back 
to fermions, the non-local OP consists of a cylindrical 
brane-type region in the bulk over which a symmetry is 
applied, and terminating operators reside on the domain 
edges. We model this general form by writing the OP 
as the long-distance limit of (OlY[j=i ^jOii)- Ol/Ob 
are possibly different operators acting near the left, right 
bosonic string (fcrmionic brane) edges and T,j is a sym- 
metry operation on a bosonic unit cell (fermionic cylin- 



drical ring). A similar form is obtained by distinguishing 
symmetry breaking bosonic phases via two-point corre- 
lations and mapping to fermions. 

Alternatively, we corroborate our conclusions by work- 
ing directly with fermions. In Appendix III, we sketch a 
derivation of fermionic selection rules. These rules deter- 
mine how the terminating operators of the fermion or- 
der parameters should be chosen, when fermion parity is 
used as the bulk symmetry, to distinguish the symmetric 
fermionic phases (Table HI)) . in analogy with the bosonic 
derivation^ The even or odd transformation rule for a 
terminating operator under P, T is sufficient and appears 
necessary to distinguish among the symmetric phases of 
any interacting model in this symmetry class. 



Phase (Bosonic variables) 


P Trans. 


T Trans. 


Example 


1. Trivial 


Even 


Even 


(Si) + 


symmetric G' — G 








2. Non-trivial 


Even 


Odd 


(S2) * 


symmetric G' — G 








3. Symmetry breaking 


Odd 


Even 


(S3) * 


G' = (T> 








4. Symmetry breaking 


Odd 


Odd 


(Si) ± 


G' = (PT) 









TABLE II. Transformation rules for Ol, Or under parity and 
translation which uniquely distinguish among the four sym- 
metric fermionic phases (N even) when fermion parity is used 
as the bulk symmetry operator. Even or odd are chosen de- 
pending on the sign of e~' M , e~ IM (Appendix III). The exam- 
ples are order parameters which are asymptotically finite in 
the listed phase and vanish in the other symmetric fermionic 
phases (primed versions S' in text also work generically) . 

Construction. Consider, as an example, an order pa- 
rameter which applies the symmetry P over many bosonic 
unit cells spanning [1,L] (large fermionic brane). Define 

i=L 

Sx = ' e i%n "y^ (19) 

The expectation (Si), taken with any choice of ground 
state, vanishes in Classes 3, 4 which have P broken in the 
bosonic variables (Appendix I). For Classes 1 and 2, P 
remains a symmetry for the bosons. Typically, applying 
a symmetry over a large domain of a symmetric state 
would yield a finite answer. 

In fact, bosonic selection rules^ can modify this con- 
clusion. The operators Ol/Or which terminate the 
bosonic string (fcrmionic brane) can be chosen to trans- 
form under symmetries in such a way as to select a 
quantum phase. The distinction between the two phases 
Classes 1 and 2 in the bosonic description are certain 
physical angles <f> = 0, 7r. O must be even under parity 
(because it is the symmetry used in the bulk of the or- 
der parameter) but transform as e 1 ^ under translation 
in order to be finite in the quantum phase labeled by 
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<fi. The OP will vanish in the other symmetric bosonic 
phase. Mapped back to the fermionic system, the brane 
termination operators should be bosonic but should be 
even or odd under translation so that the OP is finite in 
Classes 1 or 2, respectively. 

Si is, for instance, an order parameter which is finite 
in Class 1 but vanishes in Class 2 since translation in- 
variant operators terminate its bulk. To construct a 
candidate with reversed behavior, we can choose oper- 
ators such as 0(i) = *Xo(*)Xtt(*) or JXo(*)Xtt(*)j which 
are parity invariant but translation odd, to terminate 
the fermionic brane. Hence, candidate order parameters 
which uniquely give a finite value for Class 2 include 

L-l 

52^X0(1)^(1) IT e^y^xo(L)xAL) (20) 

and a similarly constructed S 2 with fermionic ends 

Xo(1)Xtt(1) and xo( L )x*{ L )- 

Finally, we utilize the fact that certain symmetries are 
broken in the bosonic variables to construct OPs that 
are finite in a symmetry breaking bosonic phase but van- 
ish elsewhere. We use two-point functions in the bosonic 
variables (UiVj) with \i — j\ — > oo; if U, V are odd under 
P but even under T, the result is finite in Class 3 but 
vanishes in the other symmetric fermionic phases. Like- 
wise, operators odd under P but even under PT yield OPs 
which can detect Class 4. Mapped back to fcrmions, a 
few such candidates are: 

L-l 

S 3 s-ixo(l) II e iwnk » (i) Xo(L) (21) 

L-l 

S 4 = -ix^l) n e^ n ^x,(L) (22) 

or similar constructions defined as S 3 , S 4 , with fermionic 
ends -ixo(l),Xo(£) and -iXir(l), Xn(L), respectively. 

Application. We apply these order parameters to 
the cylinder p + ip model. For instance, let us work 
in the t = A limit of the strong TSC with only the 
k y = DOF in the topological phase. Evaluations de- 
couple into k y = 0,7r, and fco G (0, tt) contributions: 

(S3) = {Stop)ky=0\Striv)ky=n(Striv)Y[koE(0,ir) 7^ 0- This 

OP would vanish in the other phases because the be- 
havior of the k y = 0, 7r DOF (i.e. viewed as topological 
or trivial Majorana chains) would be different. Regions 
where the order parameters take finite values are shown 
in Fig. |U 

The fcrmion parity operator used in the bulk of the 
OPs given in the previous section can be reduced to par- 
ity for just the k y — 0,7r DOF, since the remaining mo- 
menta DOF always remain trivial in the p + ip model. 
This reduction may not be applicable in general, as when 
interactions are added. 

We emphasize that this model, as in the case of the 
single Kitacv Majorana chain, has additional symmetries 
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FIG. 4. Phase diagram and non-local order parameter val- 
ues for the quasi-lD p + ip paired model (Eq. I16|l with axes 
(i,^). Phases are labeled as trivial, weak, or strong topo- 
logical superconductor (center diamond is also trivial). S/S' 
listed refer to order parameters given in text, which take on 
finite values in the phases they are associated with and vanish 
otherwise; unprimed, primed versions require either t > or 
< 0, respectively. 

beyond Z 2 x Zjv which puts constraints on the construc- 
tion of the OP. For instance, a choice oiO L / R oc (x, x) is 
different from (x, x), as we saw in the single chain. How- 
ever, as discussed previously, for models which only have 
Z 2 x Zjv symmetry, only these symmetries need to be ac- 
counted for; therefore, order parameters constructed us- 
ing the general principles described work generically. For 
instance, S3, S 3 and similar order parameters are all suit- 
able choices to distinguish Class 3 from the other sym- 
metric fermionic phases. 

Ising Condensate for Non-trivial Symmetric 
Bosonic Phase. We give a simple bosonic description 
of the non-trivial symmetric bosonic phase (Class 2) for 
the case N=2, i.e. Z 2 x Z 2 symmetry. 

Let cr, r be two Ising variables with the Z 2 symmetries 
Y\i &x(i), Yii T x(i)- Consider beginning in the partially 
ordered phase with (cr z (i)} 7^ but (r z (i)) = 0. Here, 
the domain walls of cr, created at site j by Yii<j <J x{i) 1 
and the Z 2 quanta of r, created at site j by the spin 
flip T z (j), arc both gapped. Condensing the former re- 
alizes the completely disordered phase of a and t, while 
condensing the latter results in the completely ordered 
phase. Instead, consider condensing a bound state of the 
two defects, the object p(j) = Yli<j ^M^C?)- Since the 
a domain walls have been condensed, the Ising symmetry 
of cr is restored. Moreover, though we condense r quanta, 
this does not order r since the order is non-local. Hence, 
both Z 2 symmetries are preserved. The string order pa- 
rameter for this phase is the correlation of the condensed 
objects: 

(/>(*>#)> = <r»(i) a x (k)r z (j)) (23) 

i<k<j 

The order parameter is of the general form discussed 
previously. It consists of applying one symmetry over 
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the bulk (Y[k &x(k)) and terminating with operators (t z ) 
which are even under this symmetry and odd under the 

other (YlkM*))- 

This phase is topological with protected edge states, 
which we can motivate with the following space-time pic- 
ture. In the path integral representation for this system 
with a spatial boundary, the world lines of the conden- 
sate sometimes intersect the boundary; in these cases, 
since the condensate carries a spin flip r z , the spins at 
the edges fluctuate in time. 

Based on this description, we write a Hamiltonian 
which realizes this phase. Consider starting at the critical 
point of a pair of decoupled Ising models: 

H = - ^2((J x (i) + r x (i) + a z (i)a z (i + 1) + t z (%)t z (% + 1)) 

(24) 

and adding correlations for the charge and domain wall 
bound state in order to induce condensation of the bound 
state: 

Hi = - ^2(<7 z (i)<r z (i + l)r x (i + 1) + T z (i)T z (i + l)cr x (i)) 

(25) 

The Hamiltonian H(X) = Hq + \H\ realizes the non- 
trivial topological phase for A > 1 with non-local or- 
der parameter p; this can be seen by making a dual 
transformation on one of the Z2 variables and mapping 
onto the quantum Ashkin- Teller modeh 36 ' 37 Moreover, 
Hi itself is exactly solvable and its ground state is a so- 
called cluster stated There is a four-fold degeneracy on 
a chain with sites 1 to L. On each edge, we can con- 
struct a spin-1/2 algebra with local operators; for in- 
stance, 0-3(1)7-0,(1), a x (1)tj»(2), and ct v {1)t x {1)t z {2) op- 
erate on the left edge while a x (L)r z (L), a z (L — 1)t x (L), 
and a z (L — l)a x (L)T y (L) operate on the right edge. Since 
we can map within the ground state manifold via edge 
and not bulk operators, the distinction between the de- 
generate states is topological rather than associated with 
symmetry breaking. 

V. CONCLUSION 

We have given string or brane-type non-local order pa- 
rameters to distinguish the phases of quasi-lD topolog- 
ical superconductors of spinless fermions with interac- 
tions and no symmetries. These order parameters mea- 
sure fermion parity in their bulk and are terminated by 
fcrmionic or bosonic operators at their edges; we illus- 
trated how they probe the different natures of the Majo- 
rana pairings in the topological and trivial phases. They 
would be an alternative way to detect Majorana fermions 
via quantum gas microscope measurements in cold atom 
systems. We also gave an example of an order parame- 
ter for two chains which only involves fermion parity and 
hence would be measurable using current experimental 
techniques but which detects a topological phase. 



The addition of translation to the system as a protect- 
ing symmetry distinguished among certain interesting 2D 
phases in the quasi- ID limit. We elaborated on how two 
ID Z2 invariants are distinguished by the classification 
in the case of free fermions; in particular, this allows us 
to distinguish the 2D Chern number mod 2, for instance 
the p+ip strong topological superconductor and the weak 
topological superconductor. 

We gave simple general rules which the terminating op- 
erators of a non-local order parameter should satisfy in 
order to uniquely distinguish among the fermionic sym- 
metric phases (four for N even and two for N odd) , even in 
the case of interactions. We sketched a direct fermionic 
derivation of these rules, which may be extendable to 
other symmetry classes. We illustrated the construction 
by giving a non-local order parameter for the p + ip topo- 
logical superconductor which distinguishes it from the 
weak topological superconductor or the trivial phase in 
the quasi- ID limit and which is robust to interactions. 
Attempts at extending string to brane order for coupled 
chains have been discussed in other contexts; here, we 
note that taking products of single chain string order pa- 
rameters can work because of the additional protecting 
Zjv symmetry. 

In summary, we have given uniquely identifying 
non-local order parameters for the symmetric phases of 
finitely many coupled topological superconducting spin- 
less fermion chains under interactions. We considered 
both the case of (i) no protecting symmetries and (ii) 
transverse translation symmetry. The extension of our 
ideas to incorporate time reversal symmetry, as well as 
extraction of useful 2D information, would be interesting. 
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APPENDIX 

I. Applying a Broken Symmetry on a Domain. 

Consider a quasi- ID spin system with an infinite di- 
mension indexed by j and which has a discrete bro- 
ken symmetry operator u = Y[j=-co u j- We ar g ue 
that any state |-0o) in the ground state manifold obeys 
limAr^oo^ol 11,7= -at u j\ t l J o) —> 0, with a special ordering 
of the limits. It applies even when mapped to fermions 
because it considers arbitrary ground state choices. 

Let {l^i)}^! be the broken symmetry states which 
are mapped to each other under u. We explain that 

lim A r_ i . 00 (77 i |n^L_ArWj|??fc) -s- for any i, k. If i = k, 
u creates a finite-sized domain which is orthogonal to 
the original state as the domain size increases N — > oo. 
For instance, for the quantum Ising model with a Z2 
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broken symmetry. limjv-vooCT,! IITL-Tv^jl 0, 
where | t)i I \) denote the broken symmetry states in 
the thermodynamic limit. Off-diagonal matrix elements 
i =/= k also vanish due to the order of our limiting proce- 
dures; since the thermodynamic limit precedes N — > oo 
there is always an infinite region outside the domain 
[— N, N] where the broken symmetry states are orthogo- 
nal. Practically, this means that the system size must 
be much larger than the domain over which the bro- 
ken symmetry is applied in order to yield an asymptot- 
ically vanishing value. We expect that our description 
can be formalized with matrix-product states {^qI^i 
associated with \rji) (d is the local physical dimension) 
by considering the cigenproblcm of the transfer matrices 

Eij = Ea=i K ® ( A i)*> 1 < h 3 < M, which govern the 
behavior of state overlaps. 

II. Fermionic Classification. We follow the ap- 
proach developed ir*i£ for ID fermionic and bosonic sys- 
tems. We consider the system f2 = figUfis with periodic 
boundary conditions and a unique gapped ground state. 
Let an observable be O. Consider the effective action O 
of this operator in the space spanned by the low entangle- 
ment energy (EE) Schmidt states of a subsystem fig ob- 
tained by a partition of the ground state into f2g and the 
environment He- RefJ^ observed that the action reduces 
to that of two operators Ol,O r acting locally near the 
left and right edges, respectively, of Qg, i.e. O ~ OlO r . 
That is, in this subspace spanned by low EE states, states 
are distinguished by physics near their edges (as observ- 
ables have "fractionalized" into two spatially separated 
pieces) but behave similarly in their bulk. Symmetry- 
protected phases are distinguished by the commutation 
relations obeyed by the edge operators. 

Our two Z2 x Zjv commuting symmetry generators are 
parity and translation P, T. They fractionalize as P ~ 
P L P R and f ~ T L T R . We fix P 2 = P| = P 2 = 1 and 
fN =T N =T N = 1 ( Thc proportionality constants 
depend on the following relations). Define angles fi, fi' 
with P L P R = e^P R P L , T L T R = e^'T R T L which are 
0, 7r mod 2tt since fractional pieces can be fermionic or 
bosonic. 

We claim that /n, //, along with an additional assump- 
tion that PTl = e v T^P, are sufficient to distinguish 
the quantum phases, since the other commutations fol- 
low from these. The complete operator P determines 
whether operators such as Tl are bosonic or fermionic 
(value of // for Tl) and it is natural to assume that its 
effective form docs also. Parity is in this way a more 
fundamental operator for fermionic systems relative to 
other symmetries. Other commutation relations follow, 
such as P L f = TP L e^'. 

An equation such as PlT = TP^e 1 ^ imposes a con- 
straint since P 2 = T N = 1; namely, /j,' = ir is not al- 
lowed if N is odd. Hence, we recover the same symmet- 
ric fermionic phases as in the bosonic group cohomol- 
ogy classification: (Z2) 2 for N even and Z2 for N odd. 
We additionally have a direct fermionic description of 
the phases based on effective forms. Finally, to estab- 



lish a correspondence between the bosonic and fermionic 
descriptions, we should understand when the Jordan- 
Wigner mapped versions of the fermionic symmetry op- 
erators are broken or unbroken in the bosonic variables. 
This leads us to find that parity is broken when fj, = ir, 
while translation is broken when fjf = tt and fi = 7r; this 
is summarized in Table [IJ 

III. Sketch of Fermionic Selection Rules. Wc 
sketch a proof that the terminating operators should sat- 
isfy certain transformation (selection) rules in order for 
the non-local order parameter to potentially remain fi- 
nite in one symmetric fermionic phase and vanish in 
the others; the result is Eqs. [30l We will evalu- 
ate the long-distance limit of the string or brane OP 

(Ol YijeQs R ), with Ol, O r local terminating op- 
erators and Sj an on-site symmetry. The asymptotic 
form of a non-local order parameter in a symmetric phase 
is really a two-point function of certain operators because 
symmetries reduce to acting on the edges of the domain 
over which they are applied. We use the effective forms 
for fermionic symmetries from the fermionic classifica- 
tion; though they are state dependent, we only rely on 
properties of the phases. 

We consider again a closed system f2 partitioned into a 
subsystem fig over which the symmetry £ acts and an en- 
vironment flE, on whose edges Ol,O r act. The ground 
state has Schmidt decomposition \ip) = ^ a e~ Ea \cj) a )\'q a ) 
where 4> a ,ila are for f2g,f2£, respectively. We specialize 
to the case of interest where fermion parity Pn s is ap- 
plied in the bulk. The idea ofi£ is that (<fi a \Pn s \4>a') ~ 
(4>a\Pn s ,LPris.R\ ( f ) a') (with effective forms Pn s ,L,Pn s ,R 
on f2g) for states with low entanglement energy (EE), so 
that a, a' < x with \ a cutoff. The forms Pq s ,l, Pq s ,R 
arc localized to a distance I near the edges of fig which in- 
creases with x- While the replacement by effective forms 
is approximate, it is good because states with high EE 
contribute less to evaluations of observables. Hence: 

^\0 L Pn s O R \^) « 

e- E "- E *'((t>a\Pn 3> LPn s M<t>a>)(Va\0 L R \ri a ,) 

a,a'<x 

= ]T cr E °- E «'(<f> a \(i la \0 L Pn s , L PnsMO R \(f>a>)\ria>) 

a,a'<x 

= {i>\0 L Pn Sl LPn 3 ,R0 R \i>) (26) 



where 



<X 



^a)|?7a) is a good approximation 



to ground state ip. Wc first take the thermodynamic limit 
of the closed system and then fig so that the evaluations 
at the left and right boundaries of f2g near £Ie decouple; 
the non-local order parameter reduces to an evaluation 
of local operators: 

(MOlPusOrW « (j>\0 L P ns ,L\i>){i>\Pn Sl ROR\i>) 

(27) 

We then take the limit \ J °°: s ° Pn s ,L, Pn s .R pene- 
trate further into the bulk of the (infinite) subsystem S, 
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— > and the approximation improves. 
Consider the transformation properties of just one 
edge evaluation, for instance (ip\OLPn s -L\^P)- Pn s ,L has 
known transformation rules under the symmetries, char- 
acteristic of the quantum phase; how must Ol transform 
in order to force the expression to vanish, or possibly re- 
main finite? \ip) is approximately an eigenstate of the 
effective forms of the total symmetries Pq,Tq, becoming 
exact in the above limits. Consider introducing transla- 
tion, for instance: 

{i>\0 L Pn s , L \i>) = WT^O L Pn s ,LT n \i>) (28) 
TlO L P ns , L Tn = (^ E L T aB ){J^ g P agtL T ag ) (29) 

(Note that Tq s ,Pq s are bosonic). From Appendix II, 

we have T^ g Pn s ,£Tn s = e v 'Pn s ,£. In order to keep 
the expression invariant and hence possibly finite, we 
need O l T Qb = er^' O h . Applying the same argu- 
ment with parity symmetry and using P^ Pq Si lPo s = 



e IA1 Pn s ,L implies Pq e OlPq e = er %li Oj, is needed also. 
When the end operator satisfies transformation laws dif- 
ferent from the one governed by the quantum phase 
(p,, /i') of the system, the local evaluation (^|OiPn s ,i|V') 
will vanish asymptotically. In summary, we need: 

P ] O l P = e-^0 L (30a) 
T^O L T = e- l »'0 L (30b) 

so that the non-local order parameter vanishes in the 
fermionic symmetric phases characterized by a set of an- 
gles different from (fJ.,fJ.')- These selection rules support 
the conclusions reached using bosonic selection rules and 
local order parameters for bosonic symmetry breaking. 
For instance, a non-local order parameter for Class 2 
(^, //) = (0,7r) should have Ol,Or chosen to be even 
under fcrmion parity and odd under translation, as men- 
tioned also in the main text. 
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